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A 1D2V open boundary Vlasov- Ampere code has been implemented with the aim of making a 
detailed investigation of the propagation of finite amplitude electromagnetic disturbances in an 
inhomogeneous magnetized plasma. The code is being applied to study the propagation of an exter- 
nally driven electromagnetic signal, localized at one boundary of the integration interval, through a 
given equilibrium plasma configuration with inhomogeneous plasma density and magnetic field. 

I. INTRODUCTION 

Self-consistent electromagnetic electromagnetic fields in spatially non-uniform plasmas represent one of the funda- 
mental aspects of plasma physics with several implications both in microwave and laser based experiments. Since the 
'70es extensive theoretical 0, IE an d experimental H, IE investigations have been devoted to ponderomotive 
effects, particle acceleration, wavcbreaking, resonant absorption in plasmas with density gradients. The problem is 
relevant for magnetized plasmas, as well, when Bernstein waves are excited as a consequence of mode conversion close 
to the hybrid plasma resonances. Recently, a renewed interest for electron Bernstein wave physics ha s ap peraed due 
to the possibility of implementing an attractive radiation emission diagnostic in fusion plasmas [a l& ll(T| . Then it is 
interesting to investigate the kinetic aspects of the propagation of electromagnetic as well as electrostatic fields in a 
non uniform plasma with arbitrary density and magnetic field scales, and electric field amplitudes. 

A 1D2V open boundary Vlasov-Ampere code has been implemented with the aim of achieving a more realistic 
investigation of the propagation of a finite amplitude signal in an inhomogeneous plasma, and it has been applied to 
study the propagation of an externally driven, localized charge density fluctuation in an unmagnetized plasma with 
an equilibrium spatially nonuniform density [l l| . Previous analyses of magnetized plasmas |l2l IT5 | , carried out in 
a slab geometry with periodic boundary conditions, have shown that strongly anisotropic distribution functions are 
produced both in the electron |l2j and in the ion |l3l| populations during the interaction of plasma with an externally 
applied propagating electrostatic wave. In those cases the background plasma inhomogencitics could be modelled by 
varying the ratio of the pump frequency to the electron plasma frequency |l4j . 

Here, the propagation of the spatially localized finite amplitude electromagnetic perturbation in a homogeneous 
magnetized collisionlcss plasma is prcliminarly investigated, the final aim being the study of electromagnetic fields 
propagating in a given equilibrium plasma configuration where both plasma density and magnetic field are inhomo- 
geneous. The solution of the Vlasov equation, coupled with the Ampere equation, provides the electron distribution 
function and the electrostatic field in the whole spatial range < x < L, at any time. Results are presented relevant 
to moderate amplitude electromagnetic field perturbations. 

II. THE PHYSICAL MODEL 

Let us consider a one dimensional magnetized plasma, localized in the region > x > L, homogeneous in the plane 
(U-,z), with unperturbed density uq = n e o = Zn^, and magnetic field Bo = B§e. z . The relevant dimcnsionlcss non 
relativistic Vlasov equations take the form 

f rptot TDtot~\ 

9 fa n 

where a = e,i, \x e = 1, and fa = —Zm e /mi. Normalized variables are defined as follows: t — ► tuj pe , x — ► x/d e , 

1 /2 

v — > v/c, E(B) — > eE(B)/m e cujp e . Moreover, u> pe = (47rnoe 2 /m e ) , d e = c/uj pe . As a result, in our units, the 
normalized length scale is the electron skin depth d e = 1, while the (dimensionless) electrostatic length scale, the 
Dcbye length, is equal to the normalized thermal velocity, \d — Vth.e/c. In Eq.JTJ the components of the electric field 
are the sum of two parts E l ° y = E XiV {x, t) + E^ y (x, t), the consistent field, satisfying the Maxwell equations, and the 
externally applied driving field, respectively. Moreover, the magnetic field is also the sum of the constant background 
field and of the self-consistent field, i?* ot = Bq + B z (x,t) + B z r (x,t). The self-consistent fields satisfy the Maxwell 




FIG. 1: The A;-spectra of the electric field for us — 2.1, frame (A), and us = 2.1, frame (B). 



equations: 



dE x _. dEy_ _dB 1 _ . 

at ~ Jxi at - ax Jv ' [ ' 



where the two relevant components of the current density are j x = ZriiVi X — n e V ex and j y = Zn.iViy — n e V ey . Finally, 
Poisson equation (used as a check in the code) takes the form 

dE x 

where the charge density is p = Zrii — n e . The driving electric field is modelled as 

E* r (x,t)=e 1 A; E y lr (x,t) = e 2 A; Bf(x,t) = e 2 A; A = S(x - x ) smut 

where ei f 2 are constant amplitudes, xq € [0, L] is the position where the driving field acts (usually, xq = 0), us is the 
pump frequency. The disturbance is switched on and off continuously, with a typical time scale r = 44.7. This model 
allows one to perturb the system at one boundary of the range [0, L] either by an electromagnetic (ei = 0, 62 7^ 0) or 
by an electrostatic (ei =/= 0, £2 = 0) disturbance. The consistent polarization and wave vector are then defined by the 
kinetic plasma response, which comprises the nonlinear coupling between particle motion and fields contained in the 
Vlasov Eq. |JIJ . 



III. THE RESULTS OF NUMERICAL EXPERIMENTS 



Tests aimed at reproducing the wave-plasma interaction in the low-amplitude (linear) regime have been performed 
by injecting a "pure" electromagnetic wave (i.e. E^ r = 0) at the left boundary, x = 0, with e 2 = 0.005, or by exciting 
an electrostatic perturbation (i.e. Ey T = B^ r = 0) nearby the left boundary, x ~ 0, using a normalized amplitude 
ei = 0.01. In all cases v t h. e /c = 0.14, which corresponds in dimensional units to \E XtV | 2 /(47moT e ) = 5 x 10~ 3 . In 
dimcnsionlcss units the magnetic field is equal to the electron cyclotron frequency and the value B = Q ce = 2 has been 
chosen in all simulations. Two values of the pump frequency have been considered: us = 0.95 (a) and ui = 2.1 (b). 
Note that the upper hybrid frequency ui u h = VT + B' 2 = 2.24 and the upper cutoff us co = B/2 + + B 2 /4 = 2.41. 
According to the linear theory of cyclotron waves 0| the cold branches of the dispersion relation of an extraordinary 
mode propagating perpendicularly to the magnetic field have a dimcnsionlcss wave vector k m 3.2 and k w 0.94 for 
two chosen frequencies, 2.1 and 0.95, respectively As it is seen in Fig. Q]thc fc-spectra of the electric field manifest 
a sharp maximum around the corresponding wave vectors, together with other features at higher fc's, independently 
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FIG. 2: The longitudinal, E x , and transversal, E y components of the electric field for u) = 2.1. 

of the method of excitation (electrostatic or electromagnetic). Moreover, since we are dealing with a hot plasma, 
we should expect also a component at higher k, due to the conversion of the incoming electromagnetic energy into 
electrostatic electron Bernstein waves. In this respect the most interesting spectrum is that of case (b), since the 
Cl ce < uj < u) u h, where modes with k ps 10.8 should be excited. However, in this preliminary work where we deal 
with "small" amplitude external pumps and relatively short time of propagation, this effect is not seen. In Fig. |3 
we plot the longitudinal and transversal components of the electric field, first and second frame, respectively, in the 
case of an electromagnetic external driver of amplitude = 0.005 and frequency u> = 2.1. We see that near the left 
boundary the energy is transferred to the electrostatic counterpart and then, after propagating to the right a few d e 
at x ~ 20, the amplitudes become nearly constant with E x ~ E yi despite in we are in the limit k 2 c 2 3> oj pe where 
the electromagnetic (applied) field should dominate. We think that this is due to the warm plasma response. In Fig. 
|S] we draw the isolines of the electron distribution function at t = 100 in the (x,v x ) phase space at fixed v y velocity, 
namely v v = —0.038, first frame, and v y = —0.19, second frame (B). This figure shows the rotation of the particles 
around the magnetic field corresponding vortex like structures of typical size of the order of the selected v y value. 

In the present analysis the value Vth,e/c = 0.14 has been chosen, which corresponds to an electron temperature 
of T e w lOkeV. It is well known that in order to describe correctly the collisionlcss cyclotron waves propagating 
perpendicularly to the magnetic field, the velocity dependence of the electron mass should be considered. The 
propagation of cyclotron waves is however well described already by the non relativistic theory. Moreover, for the 
frequency values which have been considered, that is uj = 0.95 and 2.1 even in the relativistic case, no appreciable 
collisionless damping is expected. 

IV. CONCLUSIONS 

In this work the Vlasov equation for electrons, coupled with the Maxwell equations, has been numerically solved in 
an open plasma slab, where the source of the electromagnetic fields is localized at one of the two spatial boundaries. 
The magnetized plasma responses to both purely electromagnetic or electrostatic excitations, at small amplitudes, 
have been studied. The field disturbance is applied in the form of a purely transverse propagating EM wave or, 
alternatively, as a purely ES field. Then, the injected EM energy goes into the longitudinal and transverse components 
of the field, almost independently on the excitation details, and the wavevectors are generated consistently during the 
wave propagation. It is seen that, even in the case of a purely EM excitation, a large fraction of the injected energy 
can go into the ES field. 

The code is particularly suitable for exploring the non linear stage of the wave-plasma interaction. The application 
of the code to large amplitude driving fields is under way. 
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FIG. 3: The electron distribution function at t = 100 in the (x,v x ) phase space with v y = —0.038, frame (A), and v y = —0.19, 
frame (B). 
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